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IX. OPTICAL PULSE PROPAGATION 

The Electromagnetic Nonlinear Schrodinger Equation: 

We begin our discussion of optical pulse propagation 35 with a derivation of the nonlinear 
Schrodinger (NLS) equation. To that end, we recall Equations [ VII- 23 ] and [ VII-23 ] 
from the early lecture set entitled Nonlinear Optics I — i.e. 

V 2 E(r,co)- V(V E(r,(o)) + (o)7c 2 e )g((o)-E(r,o))=-(i 0) 2 P (NL) (r,co) [IX-1 ] 
V (g((o) ■ E(?,co)] =- V P (NL) (r,co) [ IX-2 ] 

In this treatment we will confine our attention to wave propagation in uniform, isotropic 
optical materials - viz., glass fibers. For such materials, we can write 

P (NL) (r,Q))=8 NL (r,co)E(r,a)) [ LX-3 ] 

3 - 2 

where E NL (f ,u))= — e %^|E(r,G))| and, thus, Equation [ IX-1 ] simplifies to 

V 2 E(?,co) + (k 2 Je Q ) [e(co)+e NL (?,G))] E(?,co)= [ IX-4 ] 

where k Q = co/c- 36 

To proceed, postulate that this nonlinear Helmholtz equation can be treated by separation 
of variables methods. In particular, we are looking for a time-localized solution (a 



An excellent reference on this subject is Govind P. Agrawal's Nonlinear Fiber Optics, Academic Press (1989) 
ISBN 0-12-045140-9. 

In this simplification, we have taken V E(r,0t)) = -[e(co) -£ NL (r,G))] E(r,Q)) • V 8 NL (r,Q)) ~ 0. 
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pulse) with a relatively narrow frequency spectrum (or "group" of frequencies centered on 
a frequency 0) ctI . Thus, we assume a seoaration of variables solution 

E(f,co)= F{x,y) G(z,co-co ctr ) exp( - i |3 ctr z) [ K " 5 1 

where |3 ctI is a wave or propagation number to be associated with (0 ctI and, thus, Equation 
[ IX-4 ] becomes 

|l? + + [ t2 - p2 ] F } 5+ |0- ' 11 ,f •l r ' .' 1 5 } F = [ IX " 61 

where k= k J[e(co)+E NL (r,co)]/£ . In the linear problem (5 2 would be the "separation 

constant," but in this case we will need a bit more elaboration. Nevertheless, we shall 
assume that we can find a set of functions F(x, y) and values j3 2 that satisfy the equation 

3 2 F 3 2 F r , 2 

a7 + a7 + [ fc f 5 1 F=0 [K " 7a] 

so that 0_ /2 (3 ctr ^ + [(5 2 -(3 c 2 tr ] G=0 [ IX-7b ] 

To use perturbation theory, we first reduce Equation [ IX-7a ] to a solvable linear problem 
by writing 

k 2 = k\ [e(co)+£ NL (r,co)]/£ = [«(co)+A n(o))] 2 k\ 
~[n 2 (co) + 2n(co) An (co)] k 2 

p 2 = [(3+A(3] 2 -p 2 + 2(3A|3 [ IX-8b ] 



[ IX-8a ] 
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, e NL (r,co) 3 x£L( (D ) ifv- m 2 
where A n (©) = — — — — = — — \t,[r,(0)\ . Thus, to first order we need to solve 



2 n (u)) 8 n (u)) 



the linear equation 



3 2 F a 2 F r 2 2 2] 

3/ + 3^ + r ^"P J F=0 



[IX-9] 



which taken together with appropriate boundary conditions defines the linear eigenvalue 
problem for propagation in the medium where the functions F are the eigenfunctions and 
the values (3 are the eigenfunctions. In the previous lecture set — i.e., VIII. Guided Waves 
in Planar Structures — we found a set of eigenfunctions and eigenvalues appropriate to the 
dielectric-slab guidewave propagation 



Following an earlier discussion, we now presume that 



3 2 G 



dz 2 



az 



i.e., we take 



theslowly vary amplitude or paraxial approximation — so that Equation [ IX-7b ] 
reduces to 



3G 



where ]5 2 -p' tr -p 2 -p' tr +2(3 A(3 and, thus, 



3G 



*'2|3 ctr — +[P 2 -PL] g=o 



[IX- 10a] 



-<2|W + [P 2 -Pl] G+2(3A(3G=0 



[IX- 10b] 



For equation [ IX-7a ] to be completely satisfied in first order, we must have 
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A(3= 



_ nk^ JJ An |F {x,y)\ 2 dx dy 
P JJ |F (x,y)( dx dy 



- 3k w (3) \nl„ r.,_r., \\ 2 H| F ( X '^r ^ 



8p 



(3) I ~* / \i JJ I v '-'/I LiAtA -y 

X„,JG(z,<a-eo c „)| [K n] 



Since we are assuming that propagating pulse has a relatively narrow frequency spectrum, 
it is reasonable to use a Taylor expansion around C0 cti for (3 (co) — viz. 



PH-Pct r +(«-« c JP 1 + | ! (co-coJ 2 (3 2 + ^(co-co ctr ) 3 (3 3 + . 



[IX- 12] 



where - 



Near the dispersion minimum in glass fibers (z.e. X~ 1.3-1.6|ira ) we may, to very 
good approximation, stop with the quadratic term and write 

r 1 2 1 

p 2 - p c 2 tr = 2 f3 ctr [(3 - (3 ctr ] = 2 f3 ctr ^(co - cd ctr ) + - (co- co ctr ) (3 2 j . [IX-13] 



In this approximation, Equation [ IX- 10 ] becomes 

1 



^ + /[(co-cojP 1 + ^ 



[IX- 14a] 



If we take 



P * o 3 kl ( 3 ) , 



]\F(x,y)\ 4 dxdy 



\\F(x,y)\ dxdy 



[IX- 15] 
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Equation [ IX- 10 ] in the frequency domain reduces to 
3G 



dz 



+ i 



1 2 

(to-co c ) p\+- (co-co c ) (3 2 



g=-«y|g| 2 g 



[IX- 14b] 



which implies the following time domain equation for the pulse envelope: 
dG(z,t) dG{z,t) i d 2 G(z,t) a - 



■+P. 



dt 



~P 2 



a? 2 



+ ^G(z,?) = -/ y|g(z,?)| 2 G(z,?) [IX-16] 



The last term on the left hand side has been added to incorporate the effects of various 
possible loss mechanisms. 

Next we transform into a coordinate system which moves with the "group" — what might 
be called the surfer's coordinates of the pulse — i.e., {t,z} => {x,<;} where 
x= t— f ~Pi(z - z )and q= z- z = z-^ 1 1 . In term of these surfer's coordinates, 
Equation [ IX-16 ] becomes 37 



3%T) = q~ 
dq 2 



= - $ Gfc x) + ^ 2 - / y|G( q ,T)f % x) 



,2 



[IX-17] 



Obviously, if we omit all of the terms on the right hand side of this equation so that 
3%,x) 

— — — = 0, we would have the ideal situation wherein a pulse of any shape 
propagates forever without changing shape at a velocity v = (3 _1 -the group 



37 Since 
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velocity. In treating the less-than-ideal situation, we will initially neglect the effects of the 
loss (first) term and confine our attention to the competing effects of the dispersion 
(second) and nonlinear (third) term . We cast the lossless version of Equation [ IX- 16 ] 
into a normalized, standard form by introducing 

„(z,r) = %^, z=f = #J, r=75f 

JP Q L d Ax Ax 

where Ax is the width of the pulse and P Q is its peak power. 

Thus, we, at last, obtain the standard form of the nonlinear Schrodinger (NLS) 
equation 

du(Z,T) . , d 2 u(Z,T) 71 . ,.2 . 
i — ^ - sgn((3 2 ) -^T 1 + N 2 \u{Z,T)\ u{Z, T) = [ IX-18 ] 

where N 2 = L D /L NL = yP Ax 2 /| (3 2 1 . 

Pulse Solutions of Linear Schrodinger Equation 

If in Equation [ IX-10 ] we set G(z, i) ^> U[z,t) and neglect the loss and nonlinear terms, 
we see that U(z, t) satisfies the following differential equation: 

— V(z,t) + v- g l -U{z,t)=ib ^^z,t) [IX-19] 



. . . . . ., , 1 J 2 |3(co) 

For minimal dispersion - viz. if b = — — ■ — — 



2! did 

becomes 



= -- Equation [ IX-19 ] 



|- U(z, t) + v" 1 3- U(z, t) = [ LX-20 ] 

d z at 
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which is the basic wave equation for minimally dispersive media with the general 
solution 



where v „ = 







d(i> 


C0=C0 q _ 



U(z,t) = v(z-v g t) 
is the group velocity of the pulse. 



[ IX-21 ] 



When limitation to first order dispersion is an adequate approximation, Equation 
[ IX- 19 ] expressed in surfer's coordinates becomes 



^U(^t) + ^U(<;,t) = 



[ IX-22 ] 



where b= (3 2 /2. 



Amazingly, this pulse dispersion equation is the, so called, parabolic 
equation that we saw earlier in connection with beam propagation -- viz. 
the paraxial wave propagation equation. 

Solution of Pulse Dispersion Equation 

For convenience, we restate here Equation [ IX-22 } the first-order pulse dispersion 
equation — viz. 

b£- 2 V{q,x) + i^V{q,x) = 0. 



Let us write a Fourier transform for this modulation in terms of "surfer time" ~ i.e. 
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U(g, x) = J U(q,co) exp(/co x) dco [ LX-23a ] 

where 

U(q, co) = j- J U(q, x) exp (-z cox) Jx • [ LX-23b ] 



Thus, the pulse dispersion equation reduces to an ordinary differential equation ~ viz. 

i j- U(q,co) = b co 2 U(q,co) [ IX-24 ] 

for the Fourier transform and thus we have the simple solution 

U(^,co) = U(0,co) exp(-z b co 2 q) . [ IX-25 ] 

Thus, we see that the dispersion changes the phase of each spectral component 
of the pulse by an amount that depends on the frequency and the propagated 
distance. The general solution may be written as 

U{q, x) = J U(0,co) exp[/ (co x- b co 2 <;)] Jco [ LX-26a ] 



where 



U(0,co)= — Ju(CU- exp[-/co(?- t )] d(t- t ) ■ [ LX-26b ] 



2n _ 

Let us suppose that we have a "chirped" Gaussian at q = — i.e. 
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U(0,f- t Q )= exp 



r 



1-iC 



t-L) 



V x o J 



[ IX-27a ] 



so that 



U(0,G)): 



<i-2 ,-.2 
X n CO I 



2tc(1-iC) CXP [ 2(1- iC) J 



[ IX-27b ] 



where C > characterizes an "up-chirp" and C < a "down-chirp" pulse. 



A voltage 




A Gaussian pulse with a frequency "down-chirp" 



Inverting the transform, we see that 
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Ufc,T) = 
or rationalizing 

r 



exp<j 



(l-/C)x 2 



1 



x 2 + i2bq(l-iC)\ ^ 2[x5 + /2feq(l- iC)]j 



1 



UM= k + '2; ( i-ic)r^n]nrco 



2x; 



2 / 



l+- 



+ 



"-o y 



i+ 



2£<; v 



»o y 



xexp - 



i 




1 + C 2 ) 


X 2 


2x 2 


7 

V 




2 

+ 




2 " 



Hence, the pulse width broadening factor at a given position z is given by 
ft) 



1 + 



2£C 5 



+ 



s y 



2fr t; 



2 / .2 

1 q 



V 



+ 



^ y 



v l d y 



where L D = % 2 j2b. 



[ IX-28a ] 



[ IX-28b ] 



[ IX-29] 
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down-chirp" 



C = 



"zero-chirp" 



0.5 1 1.5 2 

RELATIVE DISTANCE (z/L D ) 



The spatial evolution of the pulse width of chirped Gaussian pulse 



(For "normal dispersion" -- i.e. 



¥1 
3co 2 



>0) 



Pulse Solutions of "Negligible Dispersion" Nonlinear Schrodinger Equation 



If group velocity dispersion can be neglected, Equation [ IX- 17 ] reduces to 



2 - 



^-5G(s.T)-/T|GfcT) -G( S .t) 



[ IX-30 ] 



and if we take G(q, x) =>£/(<;, x) yJP exp(-aq/2) we obtain 
dU{q,x) 



3? 



=- i L NL exp(-aq) | U (q,x) \ U (q,x) 



[IX-31 ] 
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where the characteristic nonlinear length is given by L NL = (y P Q ) . A solution to this 
equation is readily obtain in the form 

[ IX-32a ] 



where 



and 



U(q,x) = U(0,x) exp[/<IvM] 

<VM) = ex P(- a< ^)] [ aL NL]"' exp(-aq)| U (0,x) 
= -[q efl /h NL ]\U{0,x)\ 2 

<; eff =[l-exp(-a<;)] a" 1 



[ IV-32b ] 
[ IX-32c ] 



This interesting result show that so called "self-phase modulation" or SPM gives rise to an 
intensity dependent phase shifted or chirped pulse which remains constant 
in shape as it propagates. The instantaneous optical frequency shift is given by 



W = fa = ~LW L NlJ 



3x 



[ IV-33 ] 



Note that the pulse spectrum is "red-shifted" on the leading edge of a pulse and "blue- 
shifted" on the trailing edge of the pulse. If we suppose the initial pulse to be a super- 
Gaussian of mth-order — i.e., 



(7 m (0,T) = exp 



r l+i 



iC 



[ IV-34 ] 



- then the instantaneous optical frequency shift would be given by 



8co (x) = 



Seff 
V L NL J 



V AXy 



f x2m-l 

X 



v Ax y 



exp 



r 

V AXy 



[ IV-35 ] 
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A Brief Discussion of Solitons 

We return briefly to Equation [ IX- 18 ] — the standard form of the nonlinear 
Schrodinger (NLS) equation - i.e., 

du(Z,T) . . d 2 u(Z,T) 91 . ..2 . . 
i — ^ - s gn(P 2 ) jT 2 + l" (Z ' r) l M(Z ' r) = ° [ IX " 18 ] 
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where u (Zj) = , Z= ^ = q L5l! , r = ^ _L and JV 2 = L D /L NL = y P Q Ax 2 /| (3 2 

V^o L D Ax Ax 



In soliton analysis the most common form of NLS equation is the following: 
dU(Z,T) d 2 U(Z,T) 



dz ar 

'y Ax 2 



+ \U(Z,T)\U(Z,T) = [IX-36] 



where U(Z,T) = N u(Z,T) => G(q,x) 



P 2 



If /V > 1 the following fundamental soliton will propagate undistorted for an arbitrary 
distance: 

U(Z,T) = sech(T) exp(z Z/2) [ IX-37 ] 

If N> 2 the following second-order soliton will propagate undistorted for an 
arbitrary distance: 

4 [cosh(3 T) + 3exp (i 4 Z) cosh(r)l 

U (Z,T) = r 1 . \ , . 7^4i ex P< z / 2 [ IX " 37 ] 

v ' [cosh(4r) + 4cosh(2r) + 3cos(4Z)] FV ' ' 

]r 2 16[cosh 2 (3r) + 9cosh 2 (r) + 6cos(4 Z)cosh(3r) cosh(r)] 

\U (Z, T)\ = -j [ IX-37 ] 

[cosh(4 r) + 4 cosh(2 r) + 3cos(4 Z)\ 
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Two Soliton Collision Illustrations 



The following are time-lapsed illustration of the propagation and collision of two 
solitons ~ plot here is 

, v in 3 + 4cosh(2x-8?) + cosh(4x-64?) 
uyx, t) = 12 — — ^— ; — ^ 



[3cosh(;t-28?)+cosh(3;t -36?)] 



8 



time = -1 .00 



11 




-20 



-10 



8 



time = -1 .00 



11 




position -20 



-10 



position 



time = -0.25 


8 






6 






4 






1/ 





-10 



position 
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